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ABSTRACT 


The  objective  of  this  study  is  to  assess  quantitatively 
the  process  of  "self-development"  (Sutcliffe  and  Forsdyke, 
1950)  in  extratropical  cyclogenesis .   The  basic  approach  is 
to  determine  within  the  framework  of  the  quasi-geostrophic 
equations  the  contributions  of  thermal  advection  and  latent 
heat  release  to  changes  occurring  in  the  intensification 
rate  of  a  surface  cyclone  which  result  from  modifying  the 
field  of  vorticity  advection.   The  relevant  concepts  and 
equations  are  applied  to  a  major  cyclone  over  the  central 
United  States  and  to  a  structurally  simple  analytic  model  of 
a  cyclonic  disturbance  for  which  varying  horizontal  and 
vertical  distributions  of  latent  heat  release  are  prescribed. 

Results  of  computations  utilizing  the  analytic  model 
indicate  that  the  effect  of  latent  heat  release  in  self- 
development  can  be  comparable  with  that  of  thermal  advection. 
The  separate  and,  more  so,  the  collective  influence  of  these 
mechanisms  upon  acceleration  of  the  deepening  rate  are 
considered  synoptically  significant;  however,  the  influence 
of  translation  of  the  vorticity  maximum  aloft  with  respect 
to  the  surface  system  appears  to  dominate  changes  in  the  rate 
of  development.   The  quantitative  reliability  of  the  results 
obtained  for  the  real  data  case  are  in  doubt  because  of 
apparent  difficulties  in  the  objective  scheme  (HOBAN)  employed 
in  analysis  of  input  data. 
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I.   INTRODUCTION 

A.   OBJECTIVE 

Significant  weather  in  middle-latitudes  is  almost 
invariably  associated  with  the  development  of  extratropical 
cyclones.   Not  surprisingly,  therefore,  the  physical  and 
dynamical  mechanisms  responsible  for  cyclogenesis  have  been 
subjects  of  considerable  inquiry.   The  object  of  such  inquiry 
has  been  to  understand  more  fully  the  workings  of  extra- 
tropical  cyclones,  formulate  their  governing  laws ,  and 
accurately  predict  or  even  control  the  associated  weather. 
Since  there  are  several  excellent  sources  which  summarize 
and  synthesize  existing  concepts  of  cyclogenesis  in  middle- 
latitudes,  no  attempt  will  be  made  here  to  review  the  subject 

This  thesis  concerns  the  particular  aspect  of  extra- 
tropical  cyclogenesis  referred  to  by  Sutcliffe  and  Forsdyke 
(1950)  as  "self-development."   Specifically,  the  objective 
of  this  investigation  is  to  assess  quantitatively  the 
importance  of  the  self -development  process  and  determine  the 
relative  influence  of  the  contributing  mechanisms.   Although, 
qualitatively,  self-development  is  readily  comprehensible 
and  seemingly  important,  quantitative  assessment  of  this 
process,  to  the  best  knowledge  of  this  investigator,  has  not 
been  documented  in  the  literature. 


Of  particular  note  are  the  books  of  Petterssen  (1956) , 
Palmen  and  Newton  (19  69)  ,  and  the  thesis  of  Downey  (19  72)  . 


B.   SELF- DEVELOPMENT 

As  discussed  by  Palmen  and  Newton  (19  69) ,  self-development 
is  the  process  wherein  progressive  deformation  of  the  upper 
current  in  the  vicinity  of  a  cyclone,  as  affected  by  changes 
in  the  thermal  structure  of  the  troposphere,  results  in  a 
progressive  increase  of  the  upper-level  vorticity  advection 
and  consequent  intensified  development  of  the  surface  low. 
The  process  is  schematically  illustrated  in  Fig.  1. 
Subsequent  to  the  initiation  of  cyclogenesis ,  as  the  cyclone 
achieves  a  significant  circulation,  warm  advection  ahead  and 
cold  advection  behind  the  developing  low  center  deforms  the 
thermal  field  in  the  manner  shown  by  Fig.  1.   Deformation  of 
the  mean  isotherms  (i.e.,  thickness  lines)  through  the  layer 
considered  can  be  identified  hydrostatically  with  the  growing 
amplitude  of  the  upper-level  contour  pattern.   The  increasing 
amplitude  of  the  upper  wave  corresponds  to  an  enhanced 
variation  of  the  vorticity  (as  expressed  by  the  curvature) 
and,  therefore,  to  an  increase  of  the  vorticity  advection 
along  the  upper  current  over  the  surface  low.   This,  in  turn, 
implies  an  increase  in  the  upper-level  divergence,  which 
favors  progressively  stronger  low-level  convergence  and 
intensification  of  the  cyclone.   Note  that  at  each  instant 
of  time  thermal  advection  primarily  contributes  directly 
only  to  motion,  rather  than  to  development  of  the  sea-level 
(or  1000  mb)  system.   That  is,  at  a  given  instant,  the 
pattern  of  thermal  advection  acts  primarily  to  create  low- 
level  cyclonic  vorticity  ahead  of  the  low  and  anticyclonic 
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vorticity  behind  it,  thereby  accelerating  the  low  along  a 
path  directed  from  the  region  of  maximum  cold  to  maximum 

warm  advection.   At  that  same  time,  though,  the  thermal 

2 

advection  tends  to  enhance  the  primary  deepening  mechanism, 

namely  the  vorticity  advection  over  the  surface  low. 

If  only  horizontal  temperature  advection  were  important 
in  producing  thermal  changes  in  the  troposphere,  an  unbounded 
deformation  of  the  thermal  field  and  upper  flow  would  be 
implied.   Unbounded  increases  in  amplitude  of  the  upper  wave 
does  not  occur,  however,  because  adiabatic  cooling  and 
warming  associated  with  vertical  motions  dampen  the  influence 
of  thermal  changes  associated  with  temperature  advection.   In 
other  words,  cooling  accompanying  ascent  and  warming 
accompanying  descent  in  the  regions  of  warm  and  cold  advection, 
respectively,  tend  to  counteract  the  influence  of  horizontal 
thermal  advection.   Since  the  magnitude  of  vertical  motions 
increase  as  the  developing  circulation  intensifies,  the 
influence  of  thermal  advection  is  eventually  balanced  by  the 
opposing  influence  of  vertical  motion.   Thus,  the  process  of 
self-development  calls  into  play  a  "self-limiting"  process 
that  prevents  development  from  proceeding  without ' bound. 

In  addition  to  the  influence  of  thermal  advection  and 
vertical  motion,  thermal  changes  in  the  troposphere  can  be 
produced  by  condensational  heating.   Precipitation  and  the 


2 

Discussion  here  is  exclusive  of  the  direct  influence  of 

latent  heat  release. 
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concomitant  release  of  latent  heat  as,  for  example,  might 
occur  ahead  of  the  surface  low  center  and  in  advance  of  the 
warm  front  would  result  in  greater  amplification  of  the 
upper  ridge  downstream  from  the  low  than  would  be  the  case 
if  warm  advection  and  adiabatic  cooling  were  the  only 
processes  operative.   Self-development  would  then  proceed 
further  and  the  cyclone  achieve  greater  intensity  than  if 
precipitation  were  absent. 

C.   BASIC  APPROACH 

As  noted,  the  goal  of  this  investigation  is  to  quanti- 
tatively assess  the  process  of  self-development.   The  basic 
approach  is  to  determine,  within  the  framework  of  the  quasi- 
geostrophic  equations,  the  magnitude  and  relative  influence 
of  the  contributing  mechanisms.   More  specifically, 
determination  is  made  of  the  contributions  of  thermal 
advection  and  latent  heat  release  to  changes  occurring  in 
the  intensification  (i.e.,  deepening)  rate  of  a  surface 

cyclone  which  result  from  modifying  the  field  of  vorticity 

3 
advection.    The  relevant  concepts  and  computations  are 

applied  to  a  structurally  simple  analytic  model  of  cyclonic 

and  anticyclonic  disturbances  in  the  baroclinic  westerlies 

(after  Sanders,  1971)  and  to  a  major  cyclone  over  the  central 

United  States,  1200  GMT  10  December  1971. 


*5 

The  effect  of  vertical  motions  within  the  quasi- 
geostrophic  framework  appears  implicitly  in  tempering  the 
influence  of  thermal  advection  and  latent  heat  release,  rather 
than  explicitly  as  a  separate  contributing  mechanism. 
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II.   FORMULATION  OF  THE  PROBLEM 

A.   QUASI-GEOSTROPHIC  FRAMEWORK 

It  is  assumed  that  atmospheric  flow  is  governed  by  the 
quasi-geostrophic  vorticity  equation 


H  -  -v-v(c+f)  +  f0|f  (1) 


and  by  the  thermodynamic  equation 


|_  <|i)  -_v.v||-o(p)M-iH  ,  (2) 


where  t;  is  the  relative  vorticity;  co  is  the  vertical  velocity, 

-r^;  f   is  a  constant  value  of  the  coriolis  parameter; 

-—*-  is  the  geopotential  tendency;  R  is  the  gas  consant  for  dry 

air;  C   is  the  specific  heat  for  dry  air  at  constant  pressure; 

H  =  ■—  is  the  diabatic  heating;  and  V  =  ^—  k  x  Vcb   is  the 
dt  *  o 

geostrophic  v/ind. 

1   2 

When  the  geostrophic  relationship,  t;  =  -?—   V  §    ,  is  used 

o 

to  express  the  vorticity,  (1)  becomes 


V2  ||  =  -foWU+f)  +  f02  f   .  (3) 


r\     J. 

Combination  of  (2)  and  (3)  in  order  to  eliminate  -~   yields 

the  quasi-geostrophic  omega  equation  (see,  e.g.,  Holton,  1973), 
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2  (1) 

(aV2  +  f  2  9)W  =  f  I-  [V*V(i  V2cJ)  +  f)] 
°   8p2      °  3p       fo 

(2) 

-  y2(v*v  H> 

(3) 


p  v 


The  solution  of  this  equation  yields  values  of  to  and, 

hence,  -^—     which  enable  solution  of  (3)  for  the  field  of 
dp 

geopotential  tendency. 

The  quasi-geostrophic  equations  in  the  above  or  related 
forms  have  been  utilized  widely  (e.g.,  Sanders,  1971;  Stuart, 
1971;  Danard,  1964)  in  determination  of  fields  of  vertical 
motion  and/or  geopotential  tendency  associated  with  the 
forcing  functions  on  the  right-hand  side  of  (4) :   differential 
vorticity  advection,  term  1,  Laplacian  of  thermal  advection, 
term  2,  and  Laplacian  of  diabatic  heating,  term  3.   In  this 
investigation  emphasis  is  not  to  determine  the  direct  influence 
of  the  forcing  processes  upon  the  fields  of  vertical  motion 
or  surface  geopotential,  but  rather  to  evaluate  their  effect 
upon  sel f -development ,   That  is,  attention  is  focused  upon 
the  influence  of  the  respective  mechanisms  upon  the  field  of 
vorticity  advection  and  the  resultant  changes  that  would  be 
expected  in  the  intensification  rate  of  a  surface  cyclone. 

Application  of  (1)  to  a  level  near  the  ground  at  the 
center  of  a  cyclone  (i.e.,  at  the  vorticity  maximum  where 
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V(t;+f)  =  0)  yields  the  following  expression  for  surface 
development  in  terms  of  the  rate  of  vorticity  production  by 
low-level  convergence: 


<5C  _  f   8w 

"St    o  dp        •  (5) 


Equivalently,  one  can  write 


V  It  ~   fo  5p   •  (6) 


Here,  -rp  is  the  local  change  in  a  coordinate  system  moving 
with  the  vorticity  maximum.. 

Application  of  the  operator,  -r-r,    to  (6)  yields  an 
expression  for  the  rate  at  which  the  intensification  rate 
varies  in  time  as  a  consequence  of  changes  occurring  in  the 
field  of  vertical  motion: 


V   "St  (6t}  -  V  "6t  "   o   6p"  W  *  (7) 


Temporal  variations  in  the  field  of  vertical  motion  result 
from  changes  occurring  in  each  of  the  forcing  functions  of 
(4) .   Of  interest  in  self-development  are  variations  in  the 
vertical  motion  over  the  low  center  and,  therefore,  in  the 
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rate  of  surface  development  that  result  from  changes  in  the 

4 
field  of  vorticity  advection. 

In  order  to  evaluate  the  mechanisms  which  produce  changes 
in  the  field  of  vertical  motion  by  modifying  the  pattern  of 
vorticity  advection,  the  -ft   operator  is  applied  to  the  left- 
hand  side  and  term  (1)  of  (4) .   With  the  assumption  that  the 
stability  factor,  a,  is  invariant  in  time,  the  resulting 
expression  is 


2  (A) 

/  n2  .  r-   9   \  oo)    j-      d       r    oV    ,-, ,  1  02  .    j-\  i 

°  ^  ^    °  ^   "^     ^o"   * 

(B) 

+  foh  [v*vlt  (^  y2*)]  •         (8) 


It  can  be  seen  from  (8)  that  temporal  variations  in  the 
field  of  vertical  motion  result  from  altering  the  advecting 
velocity,  term  A,  and  the  gradient  of  relative  vorticity, 
term  B.   Modification  of  each  results  from  changes  occurring 
in  the  field  of  geopotential ,  which  can  be  obtained  from 
solution  of  (4)  and  substitution  of  the  results  into  (3). 


4 

In  a  developing  cyclone  the  Laplacian  of  thermal  advection 

over  the  low  center  is  characteristically  small  in  magnitude 
(see  e.g.,  Petterssen,  1956).   Time  variations  in  the  pattern 
of  thermal  advection  are  therefore  not  likely  to  contribute 
significantly  towards  modifying  the  rate  of  development.   On 
the  other  hand,  although  precipitation  is  typically  absent 
from  the  center  of  storms  in  the  initiation  of  cyclogenesis , 
the  Laplacian  of  diabatic  heating  over  the  surface  low  need 
not  necessarily  be  small.   Temporal  variation  in  the  amount  of 
rainfall  being  produced,  for  example,  by  a  squall  line  along 
the  cold  front  can  therefore  alter  the  rate  of  a  storm's 
development.   This  effect,  as  true  also  for  time  variation  in 
the  pattern  of  thermal  advection,  is  not  conceptually  included 
in  the  process  of  self-development  and  is  not  here  a  subject 
of  investigation. 
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The  values  of  -rr  obtained  from  (8)  on  the  basis  of  the 
ot 

computed  geopotential  tendencies  can  then  be  utilized  to 
evaluate  the  acceleration  of  the  intensification  rate, 
-T-r    (x)  t    from  (7)  .   Note  that  one  can  obtain  the  separate 
contributions  of  vorticity  advection,  thermal  advection, 
and   diabatic  heating  to  terms  A  and  B  of  (8)  and,  therefore, 
the  separate  contributions  of  each  to  -r-r-    (x)  through  the 
effects  represented  by  terms  A  and  B.   Actually,  the  influence 
of  term  A  is  not  specifically  included  in  the  concept  of 
self-development,  nor  is  the  effect  of  vorticity  advection 
on  term  B.   The  influence  of  each,  however,  is  evaluated  in 
this  study.   Figure  2  summarizes  the  mechanisms  contributing 
to  time  changes  of  the  intensification  rate  which  are  evaluated 
in  this  study. 

B.   MODEL  AND  PROCEDURES 

The  concepts  just  described  for  accessing  the  acceleration 
of  the  intensification  rate  were  applied  to  a  structurally 
simple  analytic  model  of  cyclonic  and  anticyclonic  disturbances 
in  the  baroclinic  westerlies  and  to  a  major  cyclogenesis 
over  the  central  United  States,  10  December  1971..  In  solution 
of  the  relevant  equations,  a  three-dimensional  grid  of  data 
points  is  employed.   The  vertical  grid  array  (Fig.  3)  employs 
data  at  each  100  millibar  level  from  1000  to  100  mb.   The 
horizontal  grid  measures  27  x  27  points  with  a  grid  interval 
of  160  km  at  60°  and  30°N  on  a  Lambert  Conic  Projection 
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(see  Fig.  3  for  grid  location) .   All  derivatives  are 
replaced  by  central  finite  differences,  with  the  resulting 
set  of  algebraic  equations  solved  by  the  Liebmann  relaxation 
process.   The  boundary  conditions  are  w  =  0  and  -^v-  =  0  at  all 

d  t 

lateral  and  vertical  boundaries.   The  static  stability 
factors  utilized  are  those  for  the  standard  atmosphere  as 
given  in  Table  I.   The  standard  value  of  the  coriolis  factor, 
f   is  that  for  39°N,  .92  x  10"4  sec"1. 

As  shown  in  Fig.  5,  subsequent  to  generating  the  three- 
dimensional  fields  of  geopotential  and  heating,  the  first 
step  in  the  computation  procedure  is  to  evaluate  the  contri- 
bution of  the  forcing  functions  in  (4)  to  the  field  of  vertical 
motion,  go.   The  respective  contribution  of  each  to  the 
geopotential  tendency  is  then  computed  via  (3).   Actually,  two 
sets  of  tendencies  for  each  forcing  function  are  obtained; 
the  first  being  that  associated  solely  with  the  divergence 

(-^)  term  of  (3) ,  while  the  second  is  the  total  tendency.   The 
dp 

first  set  of  tendencies  relates  to  the  developmental  influence 
of  each  forcing  function,  while  the  second  set  includes  the 
effects  of  development  and  motion. 

Once  the  field  of  geopotential  tendency  is  obtained,  the 
influence  of  each  contributing  factor  thereto  upon  -r-r   is 
evaluated  from  (8).   The  separate  contributions  of  terms  A 
and  B  are  determined.   Finally,  the  respective  influence  of 
the  relevant  physical  effects  to  temporal  variations  in  the 
rate  of  development  is  ascertained  from  (7) . 
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The  forcing  functions  of  (4)  and  ,  consequently,  all 
which  follows  therefrom  are  specified  by  the  three-dimensional 
field  of  geopotential  and  distribution  of  diabatic  heating. 
Input  for  the  real  case  study  was  height  data  at  mandatory 
data  levels  (1000,850,700,500,400,300,200,  and  100  mb) , 
converted  by  quadratic  fit  where  necessary  to  give  <f>  at  each 
100  mb  level.   Analysis  of  the  geopotential  fields  were 
accomplished  by  an  objective  technique  known  as  HOBAN 
(Inman,  1970) .   HOBAN  is  basically  a  successive  approximation 
scheme  that  starts  with  a  large  influence  radii  and  reduces 
the  radii  on  successive  passes  to  give  a  more  local  fit  on 
the  last  pass.   The  scheme  makes  use  of  height  and  wind 
within  and  outside  a  local  limited  area  grid.   The  scheme 
employs  a  distance  dependent  weighting  function  following 
Cressman  (1959).   In  an  extensive  evaluation  of  HOBAN, 
Stuart  (1974)  found  the  analyzed  height  field  and  quasi- 
geostrophic  omegas  derived  therefrom  compared  quite  favorably 
with  the  height  and  omega  fields  obtained  from  carefully 
prepared  and  consistent  subjective  analyses  that  were  used 
as  the  basis  for  judgement.   Analyses  of  the  1000  and  500  mb 
height  fields  on  1200  GMT  10  December  1971,  the  time  for 
which  computations  are  performed,  are  shown  in  Fig.  6  and 
7,  respectively. 

In  addition  to  the  real  data  case,  computations  are 
performed  for  the  structurally  simple  model  of  the  geopotential 
and  temperature  fields  of  cyclonic  and  anticyclonic 
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disturbances,  after  Sanders  (1971).   At  1000  mb ,  the 
disturbances  of  geopotential  is  given  analytically  by 


<f>(x,y,1000)  =  <j>1000  Cos  ^  (X+A)  cos  ^y  ,  (9) 


where  the  x  and  y  axes  are  directed  eastward  and  northward, 
respectively;  <f>,  oon  is  the  amplitude  of  the  pressure 
perturbation;  L  is  the  horizontal  wavelength;  and  A  is  the 
phase  lag  of  the  1000  mb  geopotential  field  relative  to  the 
temperature  field,  which  is  specified  by 


T(x,y,p)  =  TM(p)  -  (1  -  aln  i^)  (ay +  TCos  ^L  x  Cos  ^  y) 

(10) 


Here,  TM(p)  is  the  mean  temperature  at  the  pressure  level, 
p,  "a"  is  the  mean  north-south  temperature  gradient,  and 
T  is  the  amplitude  of  the  temperature  perturbation.   The 
vertical  structure  of  the  temperature  perturbation  is  governed 

by  the  factor  (1  -  aln  )  ,  which  allows  for  damping  with 

elevation  and  reversal  of  sign  at  a  pressure  level  determined 
by  the  value  of  a.   With  a  =  .722,  the  temperature  reversal 
and,  hence,  the  simulated  tropopause  occur  at  the  250  mb 
level,  as  in  middle-latitudes  in  the  real  atmosphere,  while 
the  strength  of  the  temperature  perturbation  at  500  mb  drops 
to  one-half  its  1000  mb  value.   Note  that  A  =  0  corresponds 
to  a  system  of  warm  lows  and  warm  highs.   The  1000  mb  pattern 
with  A  =  L/2  yields  cold  lows  and  warm  highs.   The  1000  mb 
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pattern  with  A  =  L/4,  representing  the  typical  intensifying 
situation  in  the  real  atmosphere,  and  with  model  parameters 
obtained  from  the  system  of  1200  GMT  10  December  1971  is 
shown  in  Fig.  8. 

At  all  levels  above  1000  mb,  the  geopotential  fields  are 
obtained  by  hydrostatic  integration;  thus, 


P      JW 

4>(x,y,p)  =  <f>(x,y,1000)  +   /   v£  (x,y,p)  dp 

1000  dp 


♦M(P)  +  ♦(xfyf1000)  -  [R  In  ^-  -    (R  |)  (In  ^-)  2  ] 


x  (ay  +  T  cos  2J2E  cos  ?IZ  )  ,  (n 


where 


1000 
4»M(p)  =  R  /    TM(p)  d  In  p  (12) 

P 


The  height  field  at  500  mb  corresponding  to  the  situation  in 
Fig.  8  is  presented  in  Fig.  9.   It  should  be  noted  that 
Fig.  8  and  Fig.  9  represent  contour  analyses  of  the  grid 
point  height  values  ascertained  from  the  analytic  expressions 
above.   Further,  all  computations  are  performed  on  the 
27  x  27  x  10  grid  in  the  same  fashion  as  for  the  real  case. 
The  main  contribution  to  diabatic  heating  is  assumed  to 
be  from  the  heat  released  by  condensation.   Other  heat  sources 
and  sinks  are  neglected.   Furthermore,  it  is  assumed  that 
condensation  occurs  as  liquid  water  which  immediately  falls 
out  as  precipitation. 
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For  the  real  data  case,  precipitation  rates  at  12  00  GMT 
10  December  1971  were  obtained  from  the  observed  two-hour 

amounts  ending  at  1300  GMT,  subjectively  averaged  over 

2 
160  km   areas  centered  at  grid  points  (Fig.  10) .   Area 

averaged  precipitation  rates  were  converted  to  the  vertically 

integrated  heating  by  multiplication  with  the  latent  heat 

factor,  L.   The  vertical  partitioning  of  the  latent  heat 

release  is  per  the  analytic  model  discussed  below,  with  the 

level  of  maximum  heating,  400  rnb ,  consistent  with  the  fact 

that  the  observed  precipitation  was  primarily  convective  in 

nature. 

In  the  analytically  specified  cyclonic  system,  the  basic 

horizontal  distribution  of  rainfall  is  modeled  as  an  ellipse 

in  which  the  maximum  precipitation  rate. P  ,  is  at  the  center 

and  values  of  P  decrease  exponentially  therefrom  (after 

Tracton,  1973) .   Thus 

2     2 

1  2     2 
A     B 
p(x,y)  =  PM  e  , 

where  A  and  B  are  scale  factors  which  specify  the  major  and 
minor  axis  (x  and  y,  respectively)  of  the  ellipse  defined 
by  P(x,y)  =  .1  PM.   The  vertical  distribution  of  heating  is 
of  the  form 

e"cz  sin  z  ,  (13) 
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where 


z  =  ( 


P  ~  Pi 
^o   1 


)7T 


Po  =  1000  mb 


p  =  200  mb  =  level  at  which  heating 
goes  to  zero   , 


and  c  specifies  the  level  of  maximum  heating. 

Computation  of  all  quantities  are  performed  for  four 
variations  of  the  basic  distribution  of  heating.   The  first, 
Case  I,  is  that  of  an  ellipse  centered  one  grid  to  the.  east 
and  one  grid  to  the  north  of  the  100  0  mb  low,  with  the  value 
of  A  and  B  such  that  the  major  and  minor  axes  are  .14  40  and 
640  km  respectively  (Fig.  11) .   The  value  of  "c"  in  this 
case  is  in  accordance  with  the  heating  maximum  in  the  vertical 
lying  at  700  mb  (Fig.  12) .   Motivation  for  the  first  case  is 
to  simulate  the  stratiform  precipitation  characteristically 
observed  ahead  of  a  low  in  advance  of  the  warm  front.   The 

maximum  intensity  of  heating  at  7  00  mb  is  explained  by  large 

dqs 
values  of  -= —   (q   =  saturation  specific  humidity)  near  the 

surface  and  large  values  of  vertical  motion  near  500  mb. 


The  correlation  between  co  and 


dp 


,  which  determines  the 


rate  of  stable  precipitation,  therefore  tends  to  be  largest 
somewhere  between  these  levels. 

The  second  case  considered,  Case  II,  is  the  same  as  the 
first,  except  that  the  level  of  maximum  heating  is  specified 
as  400  mb  (Figs.  11  and  12).   Here,  an  attempt  is  made  to 
simulate  the  effect  of  convection  embedded  in  a  warm-frontal 
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precipitation  pattern.   As  noted  by  Tracton  (1968),  warm- 
front  precipitation  is  largely  produced  by  convective 
overturning,  rather  than  by  stable  saturated  ascent.   Although 
little  is  known  about  the  detailed  distribution  of  latent 
heat  release  by  cumulus  convection,  theoretical  treatments 
of  various  investigations  (e.g.,  Kuo,  1965;  Kasahara  and 
Asai,  1967)  appear  to  agree  that  convective  heating  is 
largest  in  the  middle  to  upper  troposphere?  hence,  the  choice 
of  400  mb  for  Case  II. 

The  third  variation  of  heating  for  which  computations 
are  performed,  Case  III,  is  that  of  an  elongated  ellipse, 
1440  x  320  km,  centered  one  grid  distance  to  the  east  and 
one  grid  distance  to  the  south  of  the  1000  mb  low  and  oriented 
45°  from  the  north-south  direction  (Fig.  13) .   Motivation 
in  this  case  is  to  simulate  the  influence  of  a  squall  line 
along  or  just  ahead  of  the  cold  front.   Squall-line  precip- 
itation is,  of  course,  convectively  produced  and,  hence, 
the  level  of  maximum  heating  is  specified  to  be  at  4  00  mb. 

The  final  case  examined,  Case  IV,  represents  a  super- 
position of  the  squall  line  (Case  III)  and  the  stratiform, 
warm  frontal  precipitation  pattern  (Case  I) . 

Note  that  in  all  cases,  P„  is  specified  to  have  the 
synoptically  reasonable  value  of  .5  in  hr 
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III.   RESULTS  AND  DISCUSSION 

A.   ANALYTIC  MODEL  CASES 

Table  II  presents  the  tendencies  at  the  900  mb  low 

5 

center.    The  tendencies  are,  of  course,  a  measure  of  the 

instantaneous  rate  of  low-level  development  but  are  presented 
in  synoptically  familiar  terms  of  change  in  central  pressure 
per  hour.   Note  that  the  tendencies  associated  with 
vorticity  advection,  Xi  /  anc*  thermal  advection,  x?  '  are 
independent  of  the  diabatic  heating,  X3 '  and,  therefore, 
the  same  for  each  of  the  cases  considered. 

Exclusive  of  the  influence  of  latent  heat  release,  it  is 
clear  from  Table  II  that  vorticity  advection  is  the  dominant 
process  contributing  to  development,  i.e.,  the  active 
deepening  mechanism.   Thermal  advection,  in  accord  with  the 
qualitative  arguments  of  Petterssen  (1956) ,  produces  pressure 
falls  in  advance  and  rises  to  the  rear  of  a  cyclone,  but  has 
very  little  effect  at  the  low  center  itself  (Fig.  15) .   The 
deepening  produced  by  vorticity  advection,  -.63  mb  hr   ,  is 
consistent  with  the  value  obtained  by  Sanders  (1971) . 

The  influence  of  latent  heat  release  upon  the  deepening 
rate  varies  according  to  the  distribution  of  the  heating  with 
respect  to  the  low  center.   Values  range  from  -.55  mb  hr    in 


5 

The  900  mb  level  is  the  lowest  one  for  which  tendencies 

are  explicitly  computed.  Tendencies  at  1000  mb  required  for 
boundary  values  in  solution  of  (8)  are  obtained  by  quadratic 
interpolation . 
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the  squall  line  simulation  (Case  III)  to  -2.36  mb  hr~   in 
the  case  where  the  level  of  maximum  heating  associated  with 
warm  frontal  precipitation  is  at  700  mb  (Case  I).   Of  course, 
in  each  of  the  cases  considered  a  shift  in  the  location  of 
the  center  of  heaviest  precipitation  would  affect  the 
magnitude  of  the  computed  influence  of  heating  upon  the 
deepening  rate.   Translation  of  the  maximum  in  precipitation 
rate  to  the  low  center,  as  is  characteristic  of  the  initial 
stages  of  cyclogenesis  (Tracton,  1973),  would  result  in 
increasing  the  influence  of  condensational  heating  upon  the 
instantaneous  rate  of  development  (e.g.,  see  Fig.  16).   Note 
also  that  the  magnitude  of  the  tendencies  associated  with 
heating  are  proportional  to  the  specified  values  of  PM,  which 
in  all  the  above  cases  is  . 5  in  hr 

Values  of  the  computed  time  rate  of  change  in  intensifica- 
tion rates  at  the  900  mb  low  center  appear  in  Tables  III  and 
IV.   The  results  presented  in  Table  III  are  based  upon 
geopotential  tendencies  associated  with  only  the  divergence 
term  of  (3) ,  while  Table  IV  reflects  the  field  of  total 
tendency  produced  by  each  forcing  function.   Reference  should 
be  made  to  Fig.  2  for  the  meaning  of  terms  in  Tables  III  and 
IV.   Note  that  calculated  values  of  the  instantaneous 
acceleration  of  the  intensification  rate  have  been  applied 
to  12  hour  periods  for  ease  in  interpreting  the  results. 

The  terms  which  specifically  relate  to  the  process  of 

^2B      ^3B 
self-development  are  —tt-—   and  —^-r—    .   The  former  reflects 

dt  dt 
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the  influence  of  thermal  advection  upon  the  gradient  of 
relative  vorticity  and,  hence,  vorticity  advection  over  the 
storm  center,  while  the  latter  represents  the  effect  upon 
the  gradient  of  relative,  vorticity  resulting  from  diabatic 
heating.   Table  III  indicates  that  the  thermal  advection 
effect  produces  an  increase  in  the  deepening  rate  of  .46  mb 
hr   per  12  hr.   This  corresponds  to  an  increase  of  75  per 
cent  in  the  rate  of  development  associated  with  the  field 
of  vorticity  advection.   The  physical  mechanisms  which  are 
operative  here  are  illustrated  in  Figs.  17  to  19.   In 
response  to  the  pattern  of  thermal  advection  (Fig.  17) ,  height 
rises  occur  aloft  in  advance  of  the  surface  low,  while  height 
falls  occur  aloft  behind  the  low  position  (Fig.  18) .   The 
height  rises  and  falls  correspond  to  generation  of  anti- 
cyclonic  and  cyclonic  vorticity,  respectively.   The  enhanced 
variation  of  vorticity  and  resultant  increase  in  vorticity 
advection  along  the  upper  current  over  the  surface  low 
intensifies  the  vertical  motion  over  the  low  (Fig.  19) .   The 
consequent  increase  in  low-level  convergence  then  results  in 
an  increased  rate  of  development. 

The  influence  of  latent  heat  release  in  self-development, 
as  shown  in  Table  III,  ranges  from  -.002  to  -.50  mb  hr   per 
12  hr.   The  smallest  effect  occurs  in  the  squall  line 
simulation  (Case  III) .   This  reflects  the  fact  that  the 
height  changes  induced  aloft  by  the  heating  result  in  an 
increased  gradient  of  vorticity  across,  rather  than  along 
the  upper  flow.   Hence,  the  field  of  vorticity  advection 
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over  the  low  is  not  significantly  altered.   In  addition, 
as  can  be  seen  in  Fig.  19,  the  impact  of  the  heating  upon 
the  vertical  motion  is  principally  at  levels  above  700  mb. 
Consequently  the  effect  of  heating  upon  the  low-level 
convergence  and,  therefore,  upon  self-development  is  minimal. 

The  largest  impact  of  latent  heat  release  occurs  in  the 
700  mb  warm  front  precipitation  and  in  the  combination  700 
mb  warm  frontal  and  squall-line  precipitation  case  (Cases  I 
and  IV,  respectively) .   Note  that  the  computed  acceleration 
of  the  intensification  rate,  -.50  mb  hr~   per  12  hr  in  both 
cases,  is  comparable  in  value  to  that  ascertained  for  thermal 
advection.   Here,  the  upper  level  height  rises  induced  by  the 
heating  (Figs.  2  0  and  21)  tend  to  amplify  the  ridge  downstream 
from  the  low  center,  such  that  the  enhanced  vorticity 
advection  over  the  low  results  in  an  increased  rate  of 
intensification.   The  lesser  effect  (-.035  mb  hr   per  12  hr) 
when  the  level  of  maximum  heating  is  at  4  00  mb,  rather  than 
700  mb,  results  from  the  shift  upward  of  the  level  at  which 
heating  has  its  greatest  influence  on  the  vertical  motion. 
The  impact  of  the  heating  in  Case  IV  is  simply  not  felt  at 
low  levels. 

Recall  that  the  motivation  for  the  400  mb  warm  frontal 
precipitation  case  was  to  simulate  the  effect  of  convection 
embedded  in  large-scale  stratiform  precipitation,  which  alone 
would  have  a  level  of  maximum  heating  around  700  mb.   Actually, 
in  most  instances,  said  convection  does  not  extend  above  the 
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400  to  500  mb  level  (Tracton,  1968) ,  so  that  the  effective 

level  of  maximum  heating  would  likely  lie  somewhere  between 

3X 

700  and  400  mb .   The  total  range  of  values  of  — ~-r—  for  warm 

9t 

frontal  precipitation  is  thus  -.50  to  -.035  mb  hr~   per  12  hr. 
While  the  apparent  influence  of  convection  is  to  reduce  the 
effect  of  heating  upon  self-development  through  raising  the 
level  of  maximum  heating,  it  should  be  noted  that  the 
maximum  in  the  rainfall  rate  (to  which  all  computed  quantities 
are  proportional)  is  likely  greater  in  the  case  of  convection 
embedded  in  stratiform  precipitation  than  when  the  precip- 
itation is  purely  stratiform.   In  addition,  the  horizontal 
scale  of  heating  is  less  in  convective  situations,  so  that 
through  the  influence  of  the  Laplacian  in  the  relevant 
equations,  the  importance  of  heating  would  be  greater. 

With  reference  again  to  Case  IV,  note  that,  although  the 
combined  maximum  precipitation  rate  is  somewhat  greater  than 
that  of  Case  I  (Fig.  14  and  11,  respectively) ,  the  effect  of 
heating  on  self-development  is  exactly  the  same.   This 
reflects  the  offsetting  influence  of  an  effective  rise  in 
the  level  of  maximum  heating  and  the  shift  in  the  location 
of  the  center  of  precipitation  with  respect  to  the  low 
center. 

From  Table  III  it  can  be  seen  that  the  effect  of 
vorticity  advection  upon  the  gradient  of  vorticity  and, 

therefore,  upon  the  field  of  vorticity  advection  itself 

3Xto  _2 

(— tt— -)  ,  is  to  produce  an  increase  of  .22  mb  hr    per  12  hr  in 

0  t 
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the  rate  of  development.   As  stated  previously,  this 
mechanism  is  not  conceptually  included  in  self-development. 
The  effects  of  thermal  advection,  vorticity  advection,  and 
heating  upon  the  pattern  of  vorticity  advection  through 
their  influence  upon  the  advecting  velocity  (all  "A"  terms) 
can  be  seen  from  Table  III  to  be  negligible.   These  effects 
are  also  not  specifically  encompassed  in  self-development. 
It  can  be  seen  from  Table  IV  that  the  ^-'s  computed  on 
the  basis  of  the  total  tendency  produced  by  each  forcing 

function  are  either  negligible  in  value  or  positive  in  sign. 

3X2p 
Thus,  for  example,  the  thermal  advection  effect,   ^ 


produces  a  decrease  in  the  deepening  rate  of  .73  mb  hr 


-1 


per  12  hr,  rather  than  the  previously  discussed  increase  of 
.46  mb  hr    per  12  hr.   The  apparent  reason  for  this  is  that 
the  structure  of  the  analytical  model  is  such  that  the  maximum 
in  the  field  of  vorticity  advection  lies  ahead  of  the  surface 
low  position  (Fig.  22).   Consequently,  when  one  includes  the 
vorticity  advection  term  of  (3)  into  the  computed  tendencies, 
the  induced  translation  of  the  upper  trough  eastward  results 
in  smaller  values  of  vorticity  advection  becoming  superimposed 
over  the  surface  low.    The  resultant  reduced  rate  of  low- 
level  development  (.73- (-.46)  =  1.19  mb  hr    per  12  hr) 
is  then  sufficien  to  offset  the  opposing  influence  of  the 


c. 

The  speed  of  translation  of  the  upper  trough  is  greater 
than  that  of  the  surface  low  (Sanders,  1971). 
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the  developmental  components  of  the  forcing  functions.   It 
should  be  noted  that  the  situation  here  is  similar  to  that 
of  a  cyclone  in  the  occluding  (but  not  fully  occluded) 
stage  of  its  life  history.   In  an  earlier  stage  of  development, 
when  the  maximum  of  vorticity  aloft  lies  to  the  rear  of  the 
storm  center,  the  effect  of  translation  would  be  to  progres- 
sively increase  the  rate  of  development.   If  the  influence 
of  translation  in  such  a  case  were  of  the  same  magnitude  as 
here,  this  process  would  dominate  the  individual  contributions 
to  self-development. 

B.   CASE  OF  120  0  GMT  10  DECEMBER  19  71 

The  synoptic  case  is  that  of  a  developing  frontal  cyclone 
over  the  central  United  States  on  10  December  1971  (Fig.  6) . 
The  initial  development  of  the  storm  occurred  around  0  00  0  GMT 
10  December  over  southeastern  Oklahoma  in  association  with  an 
outbreak  of  convective  showers  in  the  vicinity  of  the  low 
center  (Fig.  23) .   By  1200  GMT,  the  convection  began  to 
spread  eastward  away  from  the  center  of  the  storm  as  the  low 
continued  to  intensify  and  move  northeastward  toward  the 
Great  Lakes.   This  association  between  the  initiation  of 
cyclogenesis  and  convective  activity,  and  the  subsequent 
dissociation  of  the  convection  from  the  low  center  has  been 
noted  by  Tracton  (1973)  to  be  characteristic  of  approximately 
one-half  the  cases  of  major  cyclogenesis  over  the  eastern 
two-thirds  of  the  United  States. 
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Figure  24  presents  a  plot  of  the  observed  change  in 
central  pressure  versus  time  for  the  period  0000  GMT  10 
December  to  00  00  GMT  11  December  1971.   The  12-hour  average 
deepening  rate  at  12  00  GMT  10  December,  the  time  for  which 
computations  are  performed,  is  1.08  mb  hr   .   The  increase 
in  the  rate  of  development  which  is  occurring  at  this  time 
is  1.67  mb  hr~   per  12  hr. 

Before  presenting  the  results  for  the  real  data  case, 
a  few  comments  are  in  order  with  regard  to  their  reliability. 
Although  the  contour  analysis  produced  by  the  HOBAN  scheme 
at  each  level  appear  synoptically  acceptable  (as,  for  example, 
Fig.  7)  and  are  generally  compatible  with  the  corresponding 
NMC  analyses,  the  vertical  consistency  between  levels  is 
somewhat  dubious.   This  is  apparent  from  the  irregularities 
in  the  derived  1000-500  mb  thickness  pattern  (Fig.  25) . 
Adjustments  to  the  HOB/\N  scheme,  which  are  beyond  the  scope 
of  this  investigation,  will  be  required  in  order  to  resolve 
the  problem.   In  addition,  Stuart  (1974)  noted  that  HOBAN 
does  well  on  the  general  pattern  and  magnitude  of  main 
centers,  but  not  so  well  on  their  precise  location.   In  the 
case  under  study  here,  this  is  manifest  in  placement  of  the 
center  of  maximum  vorticity  advection  ahead  of  the  low 
rather  than  behind  it,  as  is  true  in  the  NMC  analysis  and  is 
implied  by  the  observed  increase  in  intensification  rate. 

With  the  above  reservations,  Tables  II,  III  and  IV 
present  the  computed  tendencies  and  accelerations  of  the 
tendency  at  the  9  00  mb  low  center.   The  vorticity  advection.. 
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Table  II,  accounts  for  deepening  at  the  rate  of  -.66  mb  hr   , 
but  this  is  largely  offset  by  the  thermal  advection  term 
whose  magnitude  is  .56  mb  hr   .   HOBAN  analyses  do,  in  fact, 
show  cold  advection  at  low  levels  over  the  surface  low.   This, 
in  itself,  is  not  surprising  in  view  of  the  asymmetrical 
nature  of  the  system  in  question;  however,  there  is  little 
or  no  temperature  advection  (e.g.,  at  850  mb)  in  the 
corresponding  NMC  analysis.   Table  II  indicates  that  diabatic 
heating  contributes  negligibly  to  the  deepening  rate,  for 
the  maximum  of  the  precipitation  pattern  is  too  far  removed 
from  the  low  center  to  contribute  directly  to  its  development. 
The  total  computed  tendency  (-.12  mb  hr   )  is  clearly 
inadequate  to  account  for  the  magnitude  of  development  which 
is  observed. 

From  Table  III  one  can  see  that  thermal  advection 

contributes  -.4  3  mb  hr   per  12  hr  to  the  deepening  rate 

3^2B 
through  its  influence  upon  the  gradient  of  vorticity  (  ~   )  ; 

divergent  effect  only) .   Heating  contributes  negligibly  to 

same,  (— ^ — ),  while  the  vorticity  advection,  which  is  not 

d  t 

specifically  incorporated  in  self-development ,  produces  an 
acceleration  in  the  tendency  of  -.4  3  mb  hr    per  12  hr. 

The  combined  effect  of  the  thermal  and  vorticity  advection 
upon  changes  in  the  deepening  rate  account  for  about  one-half 
the  observed  acceleration;  however,  one  can  see  from  Table  III 
that  the  remaining  non  self-development  terms  considerably 
offset  their  influence. 
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Note  from  Table  IV  that  inclusion  of  the  vorticity 
advection  term  into  the  tendencies  required  for  computing 
acceleration  of  the  intensification  rate  renders  all  computed 
values  positive.   That  is,  all  effects  contribute  to 
decreasing  the  rate  of  development.   This  stems  from  the 
erroneous  placement  by  the  HOBAN  scheme  of  the  maximum  of 
vorticity  advection  ahead,  rather  than  behind  the  low 
center. 
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IV.   SUMMARY  AND  CONCLUSIONS 

The  objective  of  this  investigation  was  to  assess 
quantitatively  the  process  of  self -development .   Computations 
of  relevant  quantities  within  the  framework  of  the  quasi- 
geostrophic  equations  were  performed  on  a  real  data  case  and 
an  analytic  specification  of  a  cyclonic  disturbance  for  which 
varying  horizontal  and  vertical  distributions  of  latent  heat 
release  were  prescribed. 

The  quantitative  reliability  of  the  results  obtained 
for  the  real  data  case  are  dubious  because  of  apparent 
problems  in  the  objective  analysis  scheme,  HOBAN.   Clearly, 
future  effort  should  be  directed  towards  remedying,  if 
possible,  the  difficulties  in  HOBAN,  so  that  the  concepts 
and  procedures  described  herein  might  be  more  reliably 
applied  to  real  data  case  studies. 

Results  of  computations  utilizing  the  analytic  model 
indicate  that  when  only  the  divergence  (i.e.,  development) 
term  of  (3)  is  utilized  in  determination  of  the  fields  of 
geopotential  tendency  required  for  determination  of  the 
~'s,  the  effect  of  latent  heat  release  in  self-development 

dt 
^3B 

(  , ,  )  can  be  comparable  with  that  of  thermal  advection. 

ot 

Specifically,  in  both  the  700  mb  warm  frontal  precipitation 
case  (Case  I)  and  the  combined  700  mb  warm  frontal  and  squall 
line  precipitation  case  (Case  IV) ,  the  computed  time  rate 
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of  change  in  deepening  rate  due  to  condensational  heating 
was  -.50  mb  hr    per  12  hr,  versus  -.43  mb  hr    per  12  hr 
due  to  thermal  advection.   In  light  of  characteristically 
observed  accelerations  of  intensification  rate,  as,  for 
example,  the  -1.66  mb  hr    per  12  hr  for  the  10  December  1971 
case,  the  contribution  of  each  mechanisms  is  considered 
synoptically  significant,  and  their  simultaneous  occurrence 
even  more  so. 

When  the  vorticity  advection  term  of  (3)  is  included  in 
determination  of  the  fields  of  geopotential  tendency  induced 
by  the  respective  forcing  mechanisms,  the  resulting  Tr-r's 

o  t 

indicate  a  decrease,  rather  than  an  increase  in  the  rate  of 
development.   This  stems  from  the  fact  that  the  maximum  of 
vorticity  advection  in  the  model  lies  ahead,  rather  than 
behind,  the  100  0-mb  low  center,  which  is  contrary  to  the 
situation  in  the  real  atmosphere  in  association  with  a 
progressively  intensifying  cyclone.   The  magnitude  of  this 
effect,  which  represents  translation  of  the  vorticity  maximum 
aloft  with  respect  to  the  surface  system,  is  1.19  mb  hr 
per  12  hr.   Clearly,  if  the  magnitude  of  this  influence  were 
the  same  for  the  situation  where  the  maximum  of  vorticity  is 
approaching  a  surface  cyclone,  it  would  be  the  dominant 
mechanism  contributing  to  intensification  of  the  rate  of 
development. 

Note  that  the  influence  of  heating  is  dependent  upon  the 
intensity  and  three-dimensional  distribution  of  precipitation. 
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In  this  study,  the  intensity  and  distribution  are  synoptically 
and  physically  reasonable  but,  nevertheless,  somewhat 
arbitrary.   The  same  is  true  for  all  the  structural  parameters 
of  the  model,  such  as  wavelength,  intensity  of  the  temperature 
perturbation,  etc.   Further  experimentation  with  the  model 
is  therefore  in  order  to  generalize  the  results  over  a  wide 
spectrum  of  situations.   Finally,  the  concepts  and  approach 
employed  here  can  be  an  important  tool  in  diagnosing  the 
performance  of  numerical  prediction  models  in  predicting 
extratropical  cyclogenesis. 
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APPENDIX   A 


a 


944         1.178        1.2  59       2.015         3.131       4252       11.078      44.664    205.58 


Table   I.      Stability  Values,    a(m3cb    1T0N~1) 


Table    II.      Gee-potential    tendency    at    9  00    mb    low  center. 
Units:      mb   hr       (-  =    deepening,    +   =    filling) 
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Table    III.      Time    rate   of   change   of   tendency   derived   from 
divergence   term  of    (3) .      Units:      mb  hr 
(-   =    increase    in    deepening   rate) 


Table  IV.   Time  rate  of  change  of  tendency  derived  from 

divergence  and  vorticity  advection  terms  of  (3) 


Units: mb  hr"  per  12  hr  (- 
deepening  rate) 


increase  of 
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Figure  1.   Schematic  500-mb  contours  (heavy  solid  lines), 
1000-mb  contours  (thin  lines) ,  and  1000-500  mb 
thickness  (dashed) ,  illustrating  the 
"self-development"  process  during  grov/th  of  a 
cyclone.   Thermal  advection  (geostrophic)  of 
the  1000-500  mb  layer  is  inversely  proportional 
to  the  size  of  quadrilaterals  formed  by  the  1000 
contours  and  1000-500  mb  thickness  lines  (after 
Palmen  and  Nev/ton,  19  69)  . 


mb 
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Vorticity     Advection(XIA) 
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CONTRIBUTION    VIA    CHANGES    IN     THE     ADVECTING    VELOCITY 
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RELATIVE        VO  RT  IC  I  T  Y 

-i^^— -Vort  ic  ity     Advection  (Xib) 

Xtb     -* Thermal      Advection  (^2b) 

"*"     "-Heating  (X3b) 


Figure  2 .   Summary  of  contributing  factors  to  time  changes 

of  the  intensification  rate;  note  that  Xo„  ane  Xor> 

^B         oii 

are  the  only  effects  specifically  encompassed  in 
self-development . 
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0  t  Ox 

300 Z.Zm *,0m.^    .  X.ff  .  gf 


400 Z.Zm *.<?m,<«  ,   A.f^.ff 

500 Z.Zm _*i0m,w  .    X.ff.gf 

600 Z',  Zm —  0',  0m,  u   ,X,|^,| 

700 Z,Zm O.K."  .  *  .fi    .ft 
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1000 Z  ,  Zm 0,0m,W:O,     A",  ff=0 


Figure  3.   Vertical  grid,  prime  quantities  are  obtained 

by  quadratic  interpolation  from  data  at  mandatory 
levels  in  case  of  Z1  and  <$> '  ,  and  from  levels 
where  quantities  are  explicitly  evaluated  in 
case  of  X' ,  subscript  "m"  refers  to  analytic 
specification  of  initial  fields. 
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Figure  4.   Placement  of  horizontal  2  7  x  2  7  grid. 
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Figure  5.   Schematic  illustration  of  computation  procedure 


45 


±&      wo  m 


Figure  6.   HOBAN  generated  1000  mb  height  field  for 
12  00  GMT  10  December  19  71.  (NMC  frontal 
analysis  superimposed).   Units:  m. 
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Figure  7.   HOBAN  generated  500  mb  height  field  for 
1200  GMT  10  December  1971  (+  indicates 
1000  mb  low  position) .   Units:   dm. 
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Figure  10. 


Observed  rate  of  precipitation  1200  GMT 
10  December  1971.   Units:   in  hr"-1- . 
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Figure  12 . 


Vertical  distribution  of  heating;  solid  line 
represents  heating  maximum  centered  at  4  00  mb 
and  dashed  line  represents  heating  maximum 
at  700  mb. 
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Figure    19. 


Vertical  profile  of  time  rate  of  change  of 
vertical  motion,  tt-  ,  over  low  center  produced 
by  thermal  advection  (solid  line)  and  by 

diabatic  heating  of  Case  III  (dashed  line) . 

—  8        —2 
Units:   10    cb  sec 
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Figure  23.   Composite  surface-radar  summary  chart  for 
0545  GMT  10  December  1971. 
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Figure   24. 


Plot  of  Central  Pressure  (Units:  mb)  versus 
time  (Units:  hrs)  from  00  0  0  GMT  10  December 
to  0000  GMT  11  December,  1971 
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Figure  25. 


1000-500  mb  thickness  pattern  for  1200  GMT 
10  December  19  71.   Units:   dm 
(+  indicates  1000  mb  low  center) 
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